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Abstract 

We report the first observation of the decay D + — > r]'e + v e in two analyses, which combined 
provide a branching fraction of B(D + — > rj'e + v e ) = (2.16 ± 0.53 ± 0.07) x 10 -4 . We also provide 
an improved measurement of B(D + — > r]e + u e ) = (11.4 ± 0.9 ± 0.4) x 10~ 4 , provide the first form 
factor measurement, and set the improved upper limit B(D + — > cf)e + u e ) < 0.9 x 10~ 4 (90% C.L.). 
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Semileptonic decay provides an excellent laboratory for the study of both weak and 
strong interactions. Charm semileptonic decay allows determination of the parameters \V c d\ 
and \ V CS \ from the Cabibbo-Kobayashi-Maskawa (CKM) matrix pQ, and stringent testing of 
predictions for QCD contributions to the decay amplitude. A complete understanding of 
charm semileptonic decay requires study of both high-statistics and rare modes. 

The semileptonic decay D + — > 7]'e + v e has not yet been observed. Its rate relative to 
D + — > r]e + v e will provide information about rj-r]' mixing j2], as well as about the role of 
the QCD anomaly in heavy quark decays involving rf [3]. Study of these modes also probes 
the composition of the rj and rj wave functions when combined with measurements of the 
corresponding D s semileptonic decay modes jl], and can gauge the possible role of weak 
annihilation in the corresponding Z} s -meson semileptonic decays. The process D + — > <pe + v e 
is not expected to occur in the absence of mixing between the u and </>. 

The differential decay rate for D + — > 7]e + v e is given, in the limit of negligible electron 
mass, by 

rIV r 2 IT/ I 2 \rt I 3 

Ui _ ^Fl v cd\ \Plj\ 1 f (2\\2 /i\ 

where Gp is the Fermi constant, V C( i is the CKM matrix element for c — )■ d quark transitions, 
and is the r] momentum in the D meson's rest frame. The form factor f+{q 2 ) parametrizes 
the strong interaction dynamics as a function of the hadronic four-momentum transfer q 2 . 
By measuring the partial branching fraction as a function of q 2 , we probe f + (q 2 ), providing 
a test of the theoretical framework for calculation of the form factors needed to determine 
many CKM matrix elements. 

We report herein on the first observation of D + — > 7]'e + i> e and a measurement of its 
branching fraction, on an improved measurement of B(D + — > r]e + v e ) and first measurement 
of its form factor, and on an improved search for D + — > e + u e . Charge-conjugate modes are 
implied throughout this article. The results derive from two analyses of 818 pb -1 of e + e~ 
collision data collected with the CLEO-c detector [5] at the ^(3770) resonance. The data 
include ~ 2.4 x 10 6 D + D~ events. 

One analysis employs the tagging technique used in past CLEO-c studies of these [6] 
and other |8] semileptonic modes. A parent event sample is defined by reconstruction of 
either D ± meson in a specific hadronic decay mode (the tag). The fraction of these parent 
events in which the other D is reconstructed in the signal semileptonic mode determines the 
absolute semileptonic branching fraction £> SL = (-/Vta g ,SL/-^tag)(etag/etag,SL)- Here iV t ag and 
e tag are the yield and reconstruction efficiency, respectively, for the hadronic tag, and A^ ag g L 
and e tagi sL are those for the combined semileptonic decay and hadronic tag [6]. 

The six tag modes K + it~tt~, K + 7t~7t~tt , Kg7i~, Kg-K^ir , Kg7i^7i^n + , and K + K~n~ 
are selected based on the difference in energy AE = E^ — E^ cam of the D tag candidate (En) 
and the beam (-Ebeam), and on the beam-constrained mass Mbc = (-^beam/c 4 — |pz?| 2 / 'c 2 ) 1 / 2 , 
where Pd is the reconstructed momentum of the D candidate. Reference [9] summarizes 
the selection criteria and their performance for the it , K , ir°, and K$ candidates. From 
multiple candidates of the same mode and charge, we choose that with the smallest |AI£|. 
The yield of each tag mode is obtained from a fit to its M BC distribution. We find a total 
of 481223±809 D ± tags. 

We search each tagged event for an e + and an rj (77 and 7r + 7r~7r modes), rj' (7T + 7T~r] 
and p7 modes), or <j) (K + K~ mode) candidate following Ref. Candidate 7r ± e =F or K ± e T 
pairs must have an opening angle 9 > 20° to suppress 7 conversion backgrounds. For 
7/ — > p7 candidates, the ^ and 7 must have an opening angle 6 l 7r7 in the p° rest frame 
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FIG. 1. Tagged analysis U distributions in data (points) for D + — > r/'e + u e (a,b), D + — > (j)e + v e 
with <fi — > K + K~ (c), and D + — > rje + u e with r] — > 77 (d-f) and r\ — > tt + tt^tt° (g-i), in the three q 2 
intervals. The total (solid line) and background (dashed line) distributions from the fits are also 
shown. 



satisfying | cos^ 7r7 | < 0.70. Signal varies as sin 2 9 ni , while background is flat. The combined 
tag and semileptonic candidates must account for all tracks in the event. The undetected 
neutrino leads to missing energy -E m i SS = -Ebeam — Eh+ e and missing momentum p m i SS = 
-[p/ l+e + p tag (( J B b eam/c) 2 -m 2 : ,c 2 ) 1/2 ], where E h+e = E h + E e , p h+e = p^ + Pe, and p tag is the 
unit vector in the direction of the tag D~ momentum. Correctly reconstructed semileptonic 
candidates peak at zero in U = -E m i SS — c|p m i SS |, which has a 10 MeV resolution. For 
each tag mode of a given charge in an event, we allow only one semileptonic candidate. 
We take the candidate with the smallest J2x XmPOj where we sum over all reconstructed 
X G (7T°, T), rf, 4>) particles in a candidate. The pull \m{X) = (M r — Mx)/ctm, where M r 
and Mx are the reconstructed and nominal [10] masses for particle X, and the resolution 
om derives from the error matrices of the daughters of X. 

The second analysis, generic reconstruction (GR) [TT], refines techniques optimized for 
association of event-wide missing energy (-E^iss) and momentum (p^D with a neutrino [T2] . 
We apply the track and photon selection algorithms of Ref . [T2] , and impose associated event- 
level criteria to reduce background from undetected particles: the charges of the selected 
tracks must sum to zero and the number of identified e must be exactly one. We then search 
for rj (77, /t + 7t~/t°, and 37r° modes), ?/ (n + n^ri, /T 7r°^ 77 , p7, and 77 modes) candidates, using 
criteria [11] similar to those of the tagged analysis. This analysis requires | cos^Tpyl < 0.85. 

For each r]^e ± i , e candidate, the GR algorithm attempts reconstruction of a hadronic 
decay for the second D from the remaining particle content, a departure from the previous 
neutrino reconstruction measurements. Doing so both improves the E^ ss and p^L resolu- 
tions and suppresses combinatoric background. The second D reconstruction begins with 
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a closer examination of the remainder of the tracks in the event, which, again, happens 
separately for each semileptonic candidate in an event. 

From the selected tracks that are not used in the semileptonic candidate, we form two sets: 
(i) non-overlapping K s — > 7r + 7r~ candidates, and (ii) tracks consistent with originating from 
the the primary e + e~ interaction. The K s candidates must be within 12 MeV/c 2 of M K o, 
overlapping K s candidates are resolved using the best mass, and final K s candidates are 
kinematically fit with a mass constraint. A track is consistent with the primary interaction 
vertex if it is consistent with the beam envelope (within 5 cm of the origin along the beam 
direction and within 0.5 cm radially). A selected track outside of these categories is most 
likely a K s daughter whose sibling was used in the semileptonic candidate, so that candidate 
is rejected. 

To enhance photon candidate purity, we also form a set of non-overlapping 7r° — > 77 
and 7] — > 77 candidates. Overlaps are resolved based on the smallest magnitude Xm(^°) 
or Xm(v)- The algorithm's need for high efficiency dictates that we allow the broad ranges 
—25 < xm < 15 for n° candidates and — 15 < xm < 15 for r\ candidates. Unpaired showers 
with energy below 100 MeV are likely remnants from hadronic shower and are vetoed. The 
veto energy is raised to 250 MeV if any K ± candidate is found in the event. The 77 
candidates are kinematically refit with a mass constraint prior to use in the reconstruction 
of the second D and the neutrino. 

The 7T°, f], and Kg candidates, along with the remaining photons and tracks, form the 
second, non-signal (ns), D candidate with momentum p ns and energy E ns . They are further 
combined with signal e ± and 77 W candidates and compared with the total four momentum of 
the electron-positron collision to estimate E^ ss an d PmiL- The signal D momentum (p S i g ) 
and energy (E sig ) can then be reconstructed. The signal and non-signal D candidates must 
have opposite sign; the signal e ± and the non-signal D daughters must respect charge cor- 
relation assuming Cabibbo-favored decays. We require E^ ss > 50 MeV and a total vetoed- 
shower energy under 300 MeV. AE for both D candidates and E^ ss — c|p£^ ss | must be con- 
sistent with zero within mode-dependent limits of about 100 MeV. To improve resolution in 
AE sig , we take E v = c\p™i SS \- By making the further, very good, assumption that the Ip^LI 
resolution dominates the AE S [ g resolution, we can also improve Mbc- We rescale p^L^y 
a correction ( that would result in AE sig = 0: p u = CPmL w ^h ( — 1 + AE sig /(c\p^ ss \). 
Signal mode yields are determined from fits to the resulting M BC distributions. To increase 
signal sensitivity in our yield fit, we classify a high-quality (HQ) sample with the following 
properties: no unused showers, all 77 candidates with —5 < Xm(X) < 3, V G (r/,7r°), and a 
non-signal D satisfying the tagged analysis AE and Mbc criteria. Reconstruction efficien- 
cies, not including submode branching fractions, range from 2-5% overall, and 1-3% for the 
HQ subsample. 

To reduce the dominant source of background, misreconstructed decays of other more 
copious charm semileptonic modes, the GR r( candidates must satisfy Xm^D)^— x\i{.D) m ; in < 
9, where XM(-D)mm is the smallest magnitude non-signal D mass pull of all semileptonic 
candidates in an event. The additional charged and neutral D — > Xeu modes considered for 
the requirement include X = 7r ± , 7r°, i^, K* ± , K*°, K s (tt + tt~ mode), p ± , and p°. This 
requirement halves the background level with 90% signal efficiency 

Continuum backgrounds arise largely from 7 conversions or 7r° Dalitz decays in which 
one e ± lies below identification threshold. The e ± candidate is combined with each track t 
below the 200 MeV/c threshold yet with dE/dx consistent with an e T , and each e ± t T pair 
with every photon. Rejecting events with any combination satisfying m e + e - < 100 MeV/c 2 
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FIG. 2. Mbc distributions (GR analysis) for data (points) and signal (unshaded), DD (cross- 
hatch), continuum (grey), and fake (45° hatch) fit components, (a) D + — > rj'e + v e summed over 
all submodes. (b— d) D + — > r]e + u e in the indicated q 2 (GeV 2 /c 4 ) ranges, also summed over all 
submodes. 



or |m e + e - 7 — m n o\ < 50 MeV/c 2 almost completely eliminates this background. 

The 7]'e + v e yields are normalized to the K~h + 'k + yield determined using the GR tech- 
nique, but with reversal of the E^ ss requirement (-E^ss < 100 MeV) and imposition of a 
\Xm(D)\ < 3 requirement on the signal D. Other than the -E^ggrequirement, all of the 
requirements associated with the non-signal D are identical for the semileptonic modes and 
the K~ir + 7r + normalization mode. As a result, systematic effects associated with the com- 
position and reconstruction of the second D will largely cancel in the normalization ratio. 

To find q 2 = (p e + + p Ue ) 2 /c 2 , the tagged and GR analyses define the v e four momentum p v 
as (E miss , E miss p miss ) and CGpSLI, PSL), respectively. The E miss calculation in the tagged 
analysis is independent of the tag side. Using the directional information p m i SS from p m i SS 
therefore provides a more uniform calculation of q 2 across all tag modes. For the GR analysis, 
IPmLI 1S determined with better resolution than E^ ss , leading to the substitution above for 
its q 2 calculation. The D + — > r]e + v e data are divided into the q 2 ranges < q 2 < 0.5, 
0.5 < q 2 < 1.0, and q 2 > 1.0 GeV 2 /c 4 to allow study of the form factor. 

Efficiency and background determinations utilize a Monte Carlo (MC) simulation utilizing 
GEANT [13] for the detector simulation and EvtGen [H] for the physics generation. The 
analyses utilize a generic DD sample in which both D mesons decay according to the full 
model, a non-DD sample that incorporates both continuum e + e~ — >■ qq (q = u, d, or s) 
processes and radiative return production of ip(2S), and e + e~ — > r + r~ , as well as specialized 
samples for determining signal efficiency with high precision. The generic DD sample is 
equivalent to 34 times the data statistics. 

Figure [T] shows the U distributions for the tagged analysis. We observe five D + — > rf ' e r v e 
candidates: four events in the rj' — > 7t + tt~i], 77—7-77 mode, and one event in the rj' — > 7i + 7i^i], 
i] — > 7r + 7r~7r° mode. Our reconstruction efficiencies, including subsidiary branching fractions, 
are (3.26 ± 0.04)% and (0.86 ± 0.02)%, respectively, for these modes. We expect a total 
background for the combined rj' decay modes of 0.043 ± 0.026 events. Our background 
estimate is based on studies of the DD MC sample, of the generic non-DD sample, and of 
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higher statistics MC samples of background channels likely to fake the signal. The modes 
D + — > f](7T + 7i^'~f)e + i> e , D + — > ri(7T + 7T~7T°)e + i> e , and D + — > w(7r + 7r~7r°)e + ^ e with correctly 
identified tags contribute the main DD background. Using a toy simulation that folds 
Poisson statistics with statistical and systematic uncertainties, we find the probability for 
this background to fluctuate into 5 events to be 9.7 x 10 -9 , a 5.6 standard deviation (s.d.) 
significance. We find no significant signal for D + — > (f>e + v e . 

We also search for D + — > r]'e + v e with rf — > p°7 in the tagged analysis. This mode has a 
large branching fraction and detection efficiency but also a large background. No significant 
signal is observed. A 90 % C.L. upper limit is set using this decay mode: B(D + — > ?/e + z/ e ) < 
3.9 x 10~ 4 , which is consistent with the branching fractions from our observed modes. 

The rie + v e yields are determined from binned likelihood fits to the U distributions in 
each submode. The signal shape is described by a modified Crystal Ball function with two 
power-law tails [T3] that account for initial- and final-state radiation (FSR) and mismeasured 
tracks. The signal shape parameters are fixed to those determined by fits to signal MC 
samples. Background function shapes were determined by fitting the DD MC sample. 
Both normalizations float in the data fits. The main backgrounds are misreconstructed 
semileptonic decays with correctly reconstructed tags. 

The GR Mbc distributions for HQ and non-HQ samples from all rj and r\ submodes and q 2 
intervals and from K^-k + -k + are fit simultaneously with reconstructed distributions obtained 
from the MC samples for each signal mode, as well as from the generic DD and non-Dl) 
MC samples for background modeling. We employ a binned likelihood fit that incorporates 
the Barlow-Beeston methodology [16J to accommodate finite MC statistics. Simultaneous 
fitting accommodates crossfeed among all modes. The signal and DD simulations are cor- 
rected based on independent data and MC comparisons for the aspects most critical to the 
technique: the hadronic D decay model, hadronic showering in the electromagnetic calorime- 
ter, 7T° and 7] — > 77 reconstruction efficiencies, K\ energy depositions, and FSR. To probe 
the hadronic decay model, we used the GR reconstruction method with the charged and 
neutral D hadronic tags D + — > K~tt + tt + and D° — > K~tt + , respectively, in place of our 
semileptonic signal modes. We classified 108 separate decay topologies for the generically- 
reconstructed D opposite the tag. The observed rates were unfolded and efficiency-corrected, 
resulting in a decay model that, when combined with semileptonic measurements, accounts 
for 97.2% ± 2.0% of all D decays. To minimize systematic effects in this procedure, the 
rates were normalized to the unfolded D~ — > K + tt~it~ and D° — > K + tt~ to obtain branch- 
ing fraction ratios. These were then rescaled to world averages [TU] for these two modes. 
As part of this process, we also adjusted daughter spectra in the MC to reflect our data. 
The efficiency-corrected 7]'e + v e and partial r]e + u e yields float in the fit, as does the DD 
background normalization for each separate submode. Figure [2] shows excellent agreement 
between data and fit projections. 

The fit likelihood is normalized so that it would correspond to a standard \ 2 in the large 
statistics limit [17J. We find — 21n£ = 529.7 for 608 — 10 degrees of freedom, providing 
further evidence of a well-behaved fit. Fixing the r/e + z/ e yield at zero increases the — 21n£ 
by +33.6, corresponding to a statistical significance for the observed ?/e + z/ e yield of 5.8 
standard deviations. 

The statistical significance given above already incorporates both the background nor- 
malization uncertainties and finite MC statistics. Because the background normalizations 
float and the background distributions are quite flat, systematic effects must change the 
M B c shape to affect the signal yields significantly. We have used a toy MC simulation to 
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FIG. 3. Mbc distributions (GR analysis) for data (points) and signal (unshaded), DD (cross- 
hatch), continuum (grey), and fake (45° hatch) fit components for both the HQ (left) and 
non-HQ (right) subsamples in the rf — > p7 (top), r\ — > 7r + 7r _ 7/ 77 (middle), and rj — > 7r°7r°7/ 77 
(bottom) submodes. 



estimate the degradation of the significance from additive systematic effects. The toy MC 
model takes the data yields in the M BC region dominated by signal, integrated over all sub- 
modes but subdivided based on the high-quality tagging. The statistical model includes the 
independent Poisson fluctuations of the two subsamples, and the background normalization 
uncertainty that is correlated between the two subsamples. This toy model, which neglects 
some information used in the true fit, yields a statistical significance of 5.73 standard devi- 
ations, very close to our observed significance. The additive systematics are dominated by 
modeling of the K\ energy deposition, of fake charged tracks and of the momentum spectra 
of the hadronic decays. When we incorporate the additive systematic uncertainties in the 
toy MC model, taking into account the correlations between the two subsamples, we find a 
reduction in the significance that is less than 0.05 standard deviations. 

The r]'e + v e Mbc distributions for the three most influential modes are shown for both the 
HQ and non-HQ samples in Fig. |3j The data and fit are in excellent agreement across all 
these subsamples. We also examine the signal side AE, the electron momentum spectrum, 
the missing momentum spectrum and the distribution of cos Owe for both D + — > r]e + v e 
(Fig. [4| and D + — > rj'e + is e (Fig. Sj). The angle 6 We is the opening angle between the 
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FIG. 4. Comparison of data and MC components scaled by nominal GR fit results in the D + — > 
r]e + u e mode for (a) signal side AE, (b) cos 9we, (c) electron momentum spectrum, (d) missing 
momentum spectrum, (e) ir^ momentum spectrum for the non-signal D side, and (f ) tt^ momentum 
spectrum for the non-signal D side. Shown are data (points) and signal (unshaded), DD (cross- 
hatch), continuum (grey), and fake (45° hatch) fit components. 



electron and the virtual W in the W boson's rest frame, and should be distributed as 
sin 2 9we for pseudoscalar to pseudoscalar semileptonic decays such as these. The MC fit 
components are scaled according to the nominal fit results. The AE range extends outside 
of the limits imposed for the fit, and none of these distributions are used in the fit. In both 
the higher statistics D + — > r]e + v e and in the D + — > r]'e + i , e mode the scaled fit components 
and data agree very well, but would not without the signal components. The figures also 
show excellent agreement between data and the scaled fit components for the inclusive 7T 
and K ± momentum spectra from the D reconstructed against the semileptonic candidate. 
These comparisons provide strong support for our observation of D + — > r/e + z/ e . 

The systematic uncertainties in both analyses are dominated by uncertainties in the 
i] — > 77 and 7r° — > 77 detection efficiencies, with other common contributions including 
track finding efficiency, e ± , K ± and 7r =l= identification, FSR, and form- factor modeling. Effi- 
ciency and particle identification uncertainties are determined following techniques detailed 
in Ref. j9], though modified to reflect the various selection efficiencies employed by the 
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FIG. 5. Comparison of data and MC components scaled by nominal GR fit results in the D + — > 
rj'e + u e mode for (a) signal side AE, (b) cos6we, (c) electron momentum spectrum, (d) missing 
momentum spectrum, (e) momentum spectrum for the non-signal D side, and (f) momentum 
spectrum for the non-signal D side. Shown are data (points) and signal (unshaded), DD (cross- 
hatch), continuum (grey), and fake (45° hatch) fit components. 



tagged and GR analyses. The FSR and form-factor uncertainty determinations are similar 
to those in previous semileptonic analyses [6l IT2] . Other tagged contributions include un- 
certainties in iV tag , the no- additional-track requirement, and the signal U parameterization. 
The remaining GR uncertainties arise in the MC corrections described above. Many signif- 
icant uncertainties {e.g. tracking efficiency and hadronic decay model) for the GR analysis 
largely cancel in the K~-k + 'k + normalization. To account for the systematic uncertainty in 
B{D + 0e + z/ e ), we increase the upper limit by one standard deviation. 

The systematic uncertainties for both analyses are summarized in Table |TJ The r]e + u e 
quantities are signed to represent whether a given uncertainty is correlated or anti-correlated 
relative to the corresponding uncertainty for the - 0.5 GeV 2 /c 4 q 2 interval tagged analysis 
result. In forming the covariance matrix for the form factor fits for r]e + v e (see below), 
the uncertainties for a given systematic effect are treated either as fully correlated or anti- 
correlated. Treating the uncertainties from each effect as a column vector 7$, the covariance 
matrix V SJSt is then V syst = T { ® . 
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TABLE I. Systematic uncertainties (in percent) for the three D + — > r]e + v e q 2 intervals and for the 
D + — > r]'e + v e and D + — > </>e + z/ e branching fractions. The q 2 intervals are quoted in GeV 2 /c 4 . The 
r]e + u e quantities are signed to represent whether a given uncertainty is correlated or anti-correlated 
relative to the corresponding uncertainty for the 0-0.5 GeV 2 /c 4 q 2 interval tagged analysis result. 
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02 


0.02 


1.13 


1.19 


1.22 





23 


0.10 


0.60 


re 

ir —7- 77 efficiency 


0.51 





53 


0.47 


0.47 


0.36 


0.18 





75 


0.05 




77—7-77 efficiency 


3.03 


3 


07 


3.27 


2.23 


2.42 


2.87 


2 


93 


1.29 




identification efficiency 


0.66 





54 


0.43 


0.66 


0.54 


0.43 





70 


0.70 


0.80 


Simulation of FSR 


0.13 





06 


-0.15 


0.50 


0.76 


1.27 





30 


0.50 


0.30 


D + lifetime 


0.67 





67 


0.67 


0.67 


0.67 


0.67 






— 


— 


Number of D tags 


0.40 





40 


0.40 




- 







40 




0.40 


Tag fakes 


0.70 





70 


0.70 


- 


- 







70 


— 


0.70 


U fit Signal Shape 


0.37 


-0 


50 


-0.52 


- 


- 


— 






— 


— 


U fit backgrounds 


0.64 


-1 


09 


-8.25 


- 


- 


— 





22 


— 


— 


Simulation of unused tracks 


0.30 





30 


0.30 






— 





30 


— 


0.30 


Efficiency dependence on f+(q 2 ) 


1.00 


-1 


00 


-1.00 








1 


00 




3.00 


q resolution 


1.82 


-1 


90 


-0.31 
















IVlv^ blctllbllCS 


0.88 


-1 


1 n 


1 .00 








1 

1 


1U 




1 fin 

1 .uu 


K\ showering simulation 










0.00 


0.62 


1.37 






0.89 




identification efficiency 










0.05 


0.14 


0.29 






0.07 




Fake track simulation 










0.02 


0.04 


-0.08 






0.67 




Rate of unvetoed hadronic showers 










0.58 


0.05 


1.80 






0.78 




Hadronic D + decay model 










0.02 


0.04 


-0.17 






0.04 




Hadronic D + resonant substructure 










0.34 


-0.30 


0.21 






1.10 




Kg — > tt + tt~ efficiency 










0.05 


-0.12 


-0.12 






0.05 




as tt^ mis-identification 










0.02 


0.01 


-0.00 






0.00 




B(D + K-TT+7T+) 










2.20 


2.20 


2.20 






2.20 





Table IV summarizes all branching fraction and 90% confidence level (C.L.) upper limit 
results. The GR branching fractions Sgr were obtained from the measured branching ratios 
R GR = B(D+ r]'e + u e )/B(D + A^7r+7T+) using B{D + -> K-n + n + ) = (9.14±0.20)% [9]. 
The branching fractions measured using the different 77' and rj decay modes are consistent 
in both techniques. 

The tagged and GR measurements, as well as the partial f]e + i > e branching fractions within 
each measurement, are statistically and systematically correlated. To allow proper combi- 
nation of the rie + v e results, we have determined the statistical correlation matrices from 
an analysis of event overlap. Within each analysis, the statistical correlations are obtained 
from the yield fits. The statistical correlations, and combined statistical and systematic cor- 
relations (see below) are summarized in Tables |n] and III , respectively. The full correlation 
information is available from EPAPS [18] in a machine-readable format for use in fits by 
others. 
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TABLE II. Statistical correlation matrix for the three partial branching fractions for D + — > r}e + v e 
from the two analysis techniques. 





tag 


GR 




1 


-0.053 


-0.002 


0.43 








tag 


-0.053 


1 


-0.055 





0.39 







-0.002 


-0.055 


1 








0.17 




0.43 








1 


-0.043 


0.026 


GR 





0.39 





-0.043 


1 


-0.022 










0.17 


0.026 


-0.022 


1 



TABLE III. Combined statistical and systematic correlation matrix for the three partial branching 
fractions for D + — > rye + f e from the two analysis techniques.. 





tag 


GR 




1 


-0.036 


0.009 


0.439 


0.035 


0.028 


tag 


-0.036 


1 


-0.030 


0.029 


0.395 


0.018 




0.009 


-0.030 


1 


0.018 


0.014 


0.173 




0.439 


0.029 


0.018 


1 


0.030 


0.079 


GR 


0.035 


0.395 


0.014 


0.030 


1 


0.021 




0.028 


0.018 


0.173 


0.079 


0.021 


1 



To extract f+{q 2 ) for D + — > r]e + i> e , we fit the partial rates obtained from our partial 
branching fractions using td+ = 1040(7) x 10~ 15 s [10J. The fit minimizes \ 2 = A / y T V~ l A r y, 
where A7 = AT r — AT p is the vector of differences between the measured AT r and pre- 
dicted AT p partial widths, and V is the covariance matrix. We fit the two analyses both 
separately and simultaneously, taking into account statistical correlations from finite q 2 
resolution within an analysis and sample overlap between analyses. We fit first with the 
statistical covariance V = V s t a t, and then with the combined statistical and systematic 
covariance V = V stau t + V S y St . The quoted systematic uncertainties are obtained from the 
quadrature difference of uncertainties from these two fits. 

We integrate Eq. ([!]) over each q 2 interval to predict AT P , parameterizing the form factor 

TABLE IV. The branching fractions results £>tag an d Bqh from the tagged and GR analyses, 
respectively, and the branching fraction ratios i?GR relative to B{D + — > K~7r + ir + ). The errors 
are, in order, the statistical uncertainty and the systematic uncertainty. 



Mode 


Stag [10- 4 ] 


Rgk [%] 


B GR [10- 4 ] 


77V 




2.51^(0.1) 


0.237(58)(5) 


2.16(53)(7) 


<pe 




< 0.9 @ 90% confidence level (C.L.) 


rye 4 


v e 


11.1(1.3)(0.4) 


1.28(H)(4) 


11.7(1.0)(0.4) 


rye 4 


v e ,0-0.5 


6.53(94)(26) 


0.625(69)(18) 


5.71(63)(20) 


rye 4 


^,0.5-1.0 


3.08(71)(13) 


0.437(68)(13) 


3.99(62)(15) 


rye 4 


^e^l.O 


1.77(67)(16) 


0.223(52)(10) 


2.03(47)(10) 
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FIG. 6. The partial rates from the tagged (circles) and GR (squares) analyses, and the form factor 
(FF) fit (histogram). The dashed lines indicate the total uncertainty on the fit rates. 



with the standard ^-expansion parameterization [T9l [20 



We use the standard form of the outer function 0(g 2 ,t o ) and choose to to minimize the 
maximum \z\ over the physical q 2 range (see Ref. [19J). We truncate the series at k = 1 and 
allow / + (0)|V^| and the ratio of linear to constant coefficients, r\ = ai/ao to float in each 
fit. This same parameterization was used in our recent measurements of the D — > ireu and 
D — > Kev form factors [T2] • 

Figure [6] shows the combined fit, and Table V summarizes the results. For the combined 
tagged and GR fit, we find /+(0)|Kd| = 0.086 ± 0.006 ± 0.001 and r x = -1.83 ± 2.23 ± 0.28, 
with a correlation of p = 0.81. The combined fit has a x 2 — 2.5 for 4 degrees of freedom. We 
obtain the total branching fraction for the tagged and the combined analyses by integrating 
the corresponding fit result. Taking \V c d\ = 0.2256 ± 0.0010 [10J, our average value for 
|V^j|/+(0) implies /+(0) = 0.381 ± 0.027 ± 0.005. Results for other parameterizations of 
f+(q 2 ) are discussed in Appendix [XJ 

In conclusion, we have made the first observation of the decay mode D + — > r/e + z/ e and 
the first form factor determination for D + — > r/e + ^ e , as well as improving its branching 
fraction measurement. We also provide the most stringent upper limit on D + </>e + z/ e to 
date. Our combined branching fraction results are 

B(D + -> r]e + u e ) = (11.4 ± 0.9 ± 0.4) x 10" 4 , 
B(D + ->■ r]'e + u e ) = (2.16 ± 0.53 ± 0.07) x 10~ 4 , 
B(D + -»■ <j)e + u e ) < 0.9 x 10~ 4 (90% C.L.). 

These measurements are consistent with our previous results [6], which they supersede, and 
with the particle data group's upper limits |10j . They are also consistent with predictions 
from both the ISGW2 j2] and Fajfer-Kamenic [21] models. The upper limit for D + — > <pe + v e 
is about twice as restrictive as our previous limit [6]. 
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TABLE V. The D + — > rje + i , e form factor fit parameters /+(0)|V^| and n, as well as their corre- 
lation coefficient p. 



Analysis f + (0)\V cd \ n p x 2 /d.o.f. 

Tagged 0.094(9)(3) 2.17(4.50)(1.12) 0.83 0.7/(3 - 2) 
GR 0.085(6)(1) -2.89(2.24)(32) 0.81 0.0/(3 - 2) 

Combined 0.086(6)(1) -1.83(2.23)(28) 0.81 2.5/(6 - 2) 
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Appendix A: Alternate form factor parameterizations 

As our primary result for the D + — > r]e + v e form factor, we utilize the ^-expansion pa- 
rameterization [THIEn] for f + (q 2 ). This appendix provides fit results for the simple pole and 
Becirevic-Kaidalov [22J (or modified pole) parameterizations, which are commonly employed. 
The simple pole parameterization takes the form 

Utf) = r , (Al) 

pole 

while the modified pole parameterization takes the form 

h(q 2 ) = s-, (A2) 

(l-^)(l-a^ N 

Trip* ' <>> - 



where mo* is the D* mass. Either m po i e or a is fit for along with the form factor zero- 



intercept. Table VI presents the results of a combined fit to the results of the two analyses 



using these parameterizations. 
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